A NOTE ON WHY ENFORCING DISCRETE MAXIMUM 
PRINCIPLES BY A SIMPLE A POSTERIORI CUTOFF 
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Abstract. Discrete maximum principles in the approximation of partial dif- 
ferential equations are crucial for the preservation of qualitative properties of 
physical models. In this work we enforce the discrete maximum principle by 
performing a simple cutoff. We show that for many problems this a posteri- 
ori procedure even improves the approximation in the natural energy norm. 
The results apply to many different kinds of approximations including con- 
forming higher order and hp-finite elements. Moreover in the case of finite 
element approximations there is no geometrical restriction on the partition of 
the domain. 



1. Introduction 

Consider a function u : ft — > K on some bounded domain f2 e= ~R d such that 

(1.1) -Au + cusgO 

in the variational sense for some nonnegative c. Then for u + = max{0, u], we have 
the estimate 

(1.2) sup it < supu + , 

a en 

which is well-known as weak maximum principle; compare e.g. with |GT83j . 

Maximum principles usually reflect fundamental physical principles like e.g. the 
positivity of the density and it is desirable that numerical schemes respect physical 
principles. Therefore, for a conforming approximation U, computed by some nu- 
merical scheme, there arises the question if it satisfies a so called discrete maximum 
principle 

(DMP) supC/ sup£/+. 

n en 

Suppose for example that the approximation U is generated by a conform- 
ing finite element method. In this case there are plenty of results on discrete 
maximum principles; without attempting to provide a complete list, we refer to 
[CR731 IRS121 IBKK081 ILTUOl IDKS12j . All those results have in common that they 
base on piece wise affine finite element spaces and on strong geometrically restric- 
tions on the underlying partitions of the domain. To be more precise, they are 
basically restricted to non-obtuse or even acute simplicial meshes. Although there 
are some results on non-obtuse refinement in [KK111 IKK05j , it is clear that those 
restrictions introduce serious complications for the refinement and meshing of the 
domain il. The situation is even less satisfying for higher order finite elements let 
alone /ip-finitc elements. To our best knowledge for these schemes discrete max- 
imum principles are known only in a relaxed sense |Sch80j or in very restrictive 
situations; see |HM81j . 
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Based on an analysis of discrete Green's functions, Draganescu, Dupont, and 
Scott [DDS05] suggest that the discrete maximum principle for piece wise affinc 
functions may only fail in a small region close to the boundary. This motivates 
to simply cutoff those small regions with unphysical behavior. In other words, we 
define 

(1.3) U* = min{J7, sup[/+} in fl. 

en 

Obviously the modified function U* satisfies the (jDMPp independent of the poly- 
nomial degree or some underlying partition of 12. Although so far this technique 
has lacked of mathematical justification, it is actually quite popular in engineering. 

In this paper we overcome this drawback. In fact, using variational techniques 
from |DKS12[ IDLM981 IBF02] . we prove in section [2] that the modified function U* 
is even a better approximation than U. To be more precise, if u\gQ = U\gQ, then 
we have 

(1.4) \\\u-U*\\\ < |||it-f7||| 

in the norm 1 1 1 • 1 1 1 2 = § n V- 1 2 + c | • | 2 dx, induced by the reaction diffusion differential 
operator in (jl.ip . In section[3]we generalize the techniques to systems of pde's and 
the p-Laplace operator. 

We emphasize that the results are not restricted to finite element approximation 
nor to Galcrkin schemes. In fact, the results can be applied to any conforming 
approximation of u. 

2. Enforcing the Discrete Maximum Principle 

In this section we shall prove the main result (|1.3[) for functions satisfying (|1.1[) . 
To this end, for the bounded domain HcM^deN, let L 2 (Vt) and fT^fi) be the 
space of square integrable Lebesgue and Sobolev functions on 12, respectively. We 
denote by (iJ 1 (f2))* the dual of JT 1 (SI) and by -ffo(12) the subspace of functions in 
H 1 ^) with vanishing trace on the boundary 512. 

The variational formulation of (|1.1[) reads as follows: Let O^ce L 00 ^), i.e., a 
nonnegative essentially bounded function. We assume that u e if 1 (12) such that 

(2.1) Vu ■ Vw + cuvdx =: F(v) ^ for all v e Hq(Q) with v > in il. 
Jq 

With this definition we have F e (ff 1 (f2)) * and it is well known that u satisfies the 
weak maximum principle (|1.2[) ; see }GT83| . Moreover, u is the unique minimizer of 

(2.2) J(v) := ]- f | Vv\ 2 + c \v\ 2 dx - F(v) 

in u + i?g(0). In other words, u is minimal among all functions in JT 1 (17) that 
coincide with u on the boundary dQ. 

Let now U e H 1 (f)) be some approximation to u. Note that there is no restriction 
on the kind of approximation beyond of that it is conforming; we will come back 
to this issue in Remark 12.41 and the Conclusion fg] below. It follows by standard 
arguments that 

(2.3) U* := min{[/, su P C/+} = (U - sup[/+)~ + supU + e ff^fi). 

en en en 

Consequently, U* satisfies the (|DMP[) . If UigQ = u\ SS i, we have by (|1.2p that 

\u(x)-U{x)\, if U(x) < sup C/ H 

\u(x) - U*(x)\ = i x dn ^ } < \u(x) - U(x)\, 

) su P U + -u(x) 1 if U(x) > su P C/+ ' 
en en 
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for almost every ieO. Hence, it follows 

(2-4) \\u-U*\\ L2m < \\u-U\\ L2{ny 

Here || • || 2 2 ( ) '■= $ n | • | 2 da; denotes the standard norm on L 2 (f2). 

The corresponding estimate (|1 .4[) for the energy norm is less obvious. To prove 
this estimate we need some more properties of the truncated function U* . Similarly 
as in |DLM98j we have on the one hand, that 

v * i x) = { U W HU(x)< swp dQ U+, 

\snp 8Q U+ > 0, if U(x) > swp dU U+, 

and hence we obtain 

(2.5a) U* < U and \U*\^\U\ in 0. 

On the other hand we have 

[0, if U(x) >su Pan C7+, 

and consequently 

(2.5b) |V£/*| s£ |V£/| inO. 

These observations are the key properties for proving the following result. 

Theorem 2.1. Suppose the conditions of this section. In particular, let U e H 1 (fl) 
and let U* e iJ 1 (51) be the modification of U according to (|2.3[) . Then we have for 
the energy 

(2.6) and U* - U e iT x (Q). 

Proof. The second claim is a direct consequence of the definition (|2.3[) of {/* . For 
the first claim we observe from (|2.5bp that 

f |VC/*| 2 d.T < f |V[/| 2 dx. 
Jn Jn 

Furthermore, it follows from the second inequality in (|2.5ap that 

(2.7) f \U*fdx^ f |J7| 2 dx 

Jf2 JO 

and since F(v) ^ for all v e 1 (SI) , v > 0, wc obtain 

F{U) - F{U*) = F{U - U*) s£ =*> s£ --F(U"). 

Combining these estimates together with the definition (|2.2[) of the energy proves 
the theorem. □ 

Corollary 2.2. Suppose the conditions of this sections and assume further that 
U — ue Hq(CI), i.e., u = U on dfl. 

Then we have 

|||u-f7*||| < \\\u-U\\\ 
for the energy norm |||-||| 2 := ^ 2 |V-| 2 + c|-| 2 dx on H^ft) induced by (|2.1[) . 



4 



C. KREUZER 



Proof. It follows from the assumption u—Ue Hq(SY) that U,U* e u + Hq(Q). Wc 
recall that u is the unique minimizer of the energy 3 in u + Hq(Q,). Therefore, 
together with Theorem 12. 11 we have 

(2.8) 0< J(U*)- J(u) *S J(U)- J(u). 

On the other hand, for arbitrary v e u + Hq(£1) we have u — v e Hq(£V) and it 
follows from (|2.1[) that 

J(v) - J{u) = i f |Vv| 2 + c M 2 dx - - f I Vu| 2 + c |v| 2 dx + F(u - v) 
2 Jn 2 J n 

= - f |V« | 2 + c \v\ 2 dx- f • V« + cvudx + - f |Vu| 2 + c |u| 2 da; 
2 Jf2 Jn 2 J n 




Using this observation with v = U respectively v = U* in (|2.8p proves the claim. □ 

Remark 2.3. For c = t/ie maximum principle (|1.2p reads as 

sup u ^ sup u. 

n en 

Consequently, in this case, we define 

U* := xmn{U, sup[/}. 

BO. 

Since c = we do not need to have the second estimate in (|2.5a|) in order to prove 
Theorem \2.1[ All other estimates, namely the first estimate in (|2.5ap and (|2.5bl) . 
stay true for U* as defined above. Therefore, Theorem \2.1\ and Corollary \2.2\ are 
still valid in this case; compare also with the examples of Section [21 below. 

Remark 2.4. In the approximation of solutions to partial differential equations 
with finite elements one usually considers the error of the boundary values and 
the residual interior fl, separately; see e.g. )BS08] . To be more precise, let u 
satisfy (|2.1j) and let G be a discrete approximation of u on dVL. We consider an 
approximation U e 77^(0) with U = G on dtt to the weak solution uq e Hq(Q) of 
the problem 

—Aug + cuq = —Au + cu in CI, and uq = G on dfl. 

Then U satisfies Theorem \2.1\ and Corollary \2.2\ with uq instead of u, i.e., let 
U* 6 if^fi) be defined as in ([2~3]l . then 

|||u-tf*||| sS |||u g -^*||| + |||u-ug||| < |||u G -t/||| + |||«-«G||| 

The error \\\u — ug\\\ can be estimated by means of the trace u — G on 3d. Similar 
techniques can be applied in the case of a curved boundary. 



3. Extensions 

In the previous section, for the ease of presentation, we restricted ourselves to 
linear scalar valued reaction diffusion problems. However, the presented ideas can 
be generalized to more complicated problems. We shall present two examples. 
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3.1. Convex Hull Property. In this section we shall generalize the results of 
Section [5] to vector valued functions. In order to do so, we first have to generalize 
the cutoff process in (|1.3[) to higher dimensions. To this end, let K a M. m , m e N, 
be a convex set and define Hk ■ K m — * K to be the closest point projection with 
respect to the Euclidean norm | • | : K m — > R. In other words 

(3.1) U K x := argmin yeK \x - y\. 
This definition can be extended to 

U K : L 2 (fl) rn — ► L 2 (n) m setting (Il K v)(x) := U K v(x). 

It follows from the convexity of K , with elementary computations, that 11^ is 1- 
Lipschitz and hence a generalized chain rule implies Hk : H 1 (p,) m — ► i/ 1 (f2)' m 
with 

(3.2) \V(U K v)(x)\ < Lip(IIjf)|Vt;(x)| = |Vt?(x)| ie!!; 

compare with |BF021 rjD90l IAFP00] . This is the replacement for ([235) . 
Let u e H 1 ^)" 1 be such that 

| Vtt : Vv dx = for all v e Hq (0) m . 
Jn 

It is well known that u satisfies the convex hull property 

u(VL) a convhullu,((7Sl) 
(see |BF02| ') and that it is the unique minimizer of the energy 

J(v) := \ \\Vv\ 2 dx mu + H^(n) m . 
Jn 2 

Let U e 1 (Sl) m be some approximation of u and define 

U* := U K U, with K := convhull U(dCl). 

Consequently U* satisfies the discrete convex hull property 

U* (fi) c convhull 17* (dO) 

and it follows from (|3.2|) similar as in the proof of Theorem 12.11 that U* — U e 
H^(fl) m and 

(3.3) J(U*) = l - f \VU*\ 2 &x^ 1 - [ \VU\ 2 dx = J(U). 

2 Jn 2 J n 

Moreover, as in Corollary 12.21 we obtain for |||-||| 2 = § n |V-| 2 da; that 

(3.4) |||u - U*||| < |||u- U\\\ if u-UeH^(Q) m . 

Remark 3.1. The results of this section can be generalized to reaction diffusion 
problems. To this end we observe that, in order to prove (|3.3[) . we need a replace- 
ment of (|2.5aj) for vector-valued functions. Note that for reaction diffusion problems 
the convex hull property reads as 

u(O) c convhull ({0} u u(8fl}). 

Therefore, defining K := convhull({0} u U(dCl)) we obtain \U (x)\ ^ |l/(x)|, 
x e Q. This is the required replacement of (|2.5ap . 
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3.2. Nonlinear Problems. Recalling the proof of Theorem 12.11 and (|2.5p . we ob- 
serve that formally (|2.6p holds for energies 

J{v) = f T(x,\v\,\Vv\)dx- F(v) 
Jn 

with J:!]xIxl^I being monotone in its second and third variable. In 
particular, Theorem 12 . 1 1 can be applied to many nonlinear problems. 

As an example, we consider the nonlinear p-Laplacc problem. To this end, for 
fixed p £ (1, oo), i + i = 1, let W 1,p (fi) be the space of p integrable functions on SI 
with p-intcgrable weak derivatives. Let u e W 1,p (£l) such that 

-div|Vu| p ~ 2 Vu =: F < in SI; 

in the distributional sense. Then u is the unique minimizer of the energy 

J{v):= \ -\Vv\ p dx- F{v) in u + W 1,p (ft). 

Here Wo' p (f2) is the subspace of functions in W 1,P (Q) with vanishing trace on dfl. 
Moreover, it is well known that u satisfies 

sup u ^ sup u; 

n en 

see e.g. [DLM98] . Let U e W 1,V (Q) be some approximation to u and let 

U* : = min { U, sup U}. 

en 

This implies (|2.5bp . Note that the differential operator contains no reactive term 
and thus we do not require (|2.5aj) to conclude, similar as in the proof of Theorem 
|2~T1 that 

J(U*)^J(U) and thus J(U*) - J{u) < J(U) - J{u). 

In order to prove a result analog to Corollary 12.21 it remains to correlate the 
energy difference to a reasonable measure of distance. This can be done using 
the so-called quasi- norm introduced by Barrett and Liu in |BL93j . It follows from 
|Krel2[ Lemma 13] that there exist constants C, c > such that 

c(j(v) - J{u)) ^\\\v- u\\\\ Vu) := J (M + \v\Y' 2 \u - v\ 2 dx < C(j(v) - J{u)) 

for all v e u + W ' P (Q). Hence, with v = U respectively v = U*, it follows from 
above, that 

V~c\\\u-U*\\\ {Vu) ^VC\\\u-U\\\ {Vu) . 
4. Conclusion 

Maximum principles often reflect physical behavior of solutions and it is there- 
fore desirable that numerical approximations satisfy a maximum principle as well. 
In this paper we enforce the maximum principle by performing a simple cutoff to 
the approximation. We show that for many problems this truncation even im- 
proves the approximation. Hence, all error estimates and convergence results for 
the approximation can directly be applied to the error of the truncated function. 

We emphasize that we do not specify the sort of approximation beyond that it is 
conforming. Therefore, among others, all presented results apply to all conforming 
finite element approximation including high order elements and conforming hp- 
methods. There is no restriction on the underlying partition of SI like non obtuseness 
of the triangulation and the results even apply to partitions involving complicated 
element geometries. 
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